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In the present work, we have considered N-dimensional Einstein field equations in which 
4-dimensional space-time is described by a FRW metric and that of the extra d-dimensions 
by an Euclidean metric. Considering the universe filled with tachyonic field we have recon- 
structed the potential V((f>) corresponding to the field reconstructions in an anisotropic 
O 

universe under Emergent-Power law, Emergent-Intermediate, Emergent-Logamediate and 



5-H 



> 



Logamediate-Intermediate scenarios. The statefinder parameters have been investigated in 
all of the said scenarios. 



INTRODUCTION 



m 

The two independent observational signals on distant Type la supernovae (SNIa) in 1998 have 
revealed the speeding up of our universe. This acceleration implies that if the theory of Einstein's 
gravity is reliable on cosmological scales, then our universe is dominated by a mysterious form 
of matter. This unknown component possesses some remarkable features, for instance it is not 

a ■ 

clustered on large length scales and its pressure must be negative in order to be able to drive the 
current acceleration of the universe [l(. This matter content is dubbed as "dark energy" (DE) that 
occupies about 70% of today's universe. Reviews on DE include 
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A chronological 



history of the research on DE is discussed in [J] . The basic characteristic of DE is that its equation 
of state (EOS) parameter w = p/p, where p is the energy density and p is the pressure that has a 
negative value. The simplest candidate of dark energy is a tiny positive cosmological constant [2] 
corresponds to a fluid with a constant equation of state w = —1. However, as is well known, it is 
plagued by the so-called "cosmological constant problem" and "coincidence problem" [2(]. Other 



'Electronic address: chayanranjit@gmail.com 

^Electronic address: surajit'2008@yahoo.co.in, surajit.chattopadhyay@redifTmail.com 



dark energy models include quintessence [5J], phantom [6], quintom [1(, Chaplygin gas [7J], tachyon 
a\, hessence [9|] etc. 

The weakness of the gravitational force has been successfully explained by postulating the exis- 
tence of extra dimensions [10l |. Today, there is a large variety of promising models and theories 
which suggest the existence of more than three spatial dimensions. Most notably, string theory 
suggests the existence of seven additional spatial dimensions [111 ]. As our space-time is explicitly 
four dimensional in nature so the 'hidden' dimensions must be related to the dark matter and dark 
energy which are also 'invisible' in nature. Model of higher dimension was proposed by Kaluza and 
Klein |12j, |l3j who tried to unify gravity with electromagnetic interaction by introducing an extra 
dimension which is basically an extension of Einstein general relativity in 5D. The activities of ex- 
tra dimensions also verified from the STM theory [14| proposed recently by Wesson et al. [15|] . The 
importance of extra dimension in cosmology has been discussed by many authors. In some analysis 
the metric is assumed to be uniform in extra dimension [18|, [l9|]. For the metric having non-trivial 



dependence the extra dimension, the analysis on inflation may be altered 



20, 



211 ]. Panigrahi and 



Chatterjee 22| have found that the inflationary scenario is possible for inhomogeneous extra di- 



mensional model and for the homogeneous case an initially decelerating universe starts accelerating 



undergoing a flip. Ibdnez and Verdaguer 



231 ] have obtained a set of solutions of Einstein's equation 



in an N-dimensional vacuum model and also homogeneous solutions with expanding 3-dimensional 
isotropic space. The 4D perfect fluid solutions in flat universe model have been obtained by Krori 
et al. [24|, Gleiser and Diaz [251 ] . using a higher dimensional anisotropic cosmology (Bianchi-I) 
which are compatible with contraction of all the extra dimensions. In [261 ] Paul has considered the 
theories of imperfect fluid such as Eckart, EIT, TIS and FIS to obtain the cosmological solutions 
for flat FRW with extra dimensions by Kasner type Euclidean metric and present an analysis of a 
n-dimensional vacuum Einstein's field equations. In [27] Gorbunov and Sibiryanov have proposed 
a cosm olog ical model of self accelerated brane universe with warped extra dimension. In Peng 



et al. 



281 ] extended the direct quantum approach to the FRW cosmology from 4D to 5D and 



29] 



obtained a Hamiltonian formulation for a wave like 5D FRW cosmology. In Panigrahi et al. 
have shown a scenario in-homogeneous 5D space time which behave a decelerating expansion in 
the early epoch along with an accelerating situation at the present line without introducing any 
external quintessence-like scalar field in the presence of extra dimension. 

In the present work, we have considered N-dimensional Einstein field equations in which 4- 
dimensional space-time is described by a FRW metric and that of the extra d-dimensions by 
an Euclidean metric. Considering the universe filled with tachyonic field we have reconstructed 



the potential V(cf>) corresponding to the field reconstructions in an anisotropic universe. A rolling 
tachyon has an interesting equation of state whose parameter smoothly interpolates between —1 
and 30|] . This has led to several attempts to construct viable cosmo 



tachyon as a suitable candidate for the inflaton at high energy 

flat FRW background the energy density and pressure of tachyon are given by [2 ] 
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where V and <j> denote the potential and tachyonic field respectively. The tachyonic matter might 
provide an explanation for inflation at the early epochs and could contribute to some new form 
of cosmological dark matter at late times 



321 ] . In isotropic FRW universe, the tachyonic field 
was studied by [33j], where it was observed by studying the statefinder diagnostics [34f] that the 
tachyonic interpolates between dust and ACDM stages of the universe. In the present paper we 
are going to consider the scale factors a{t) and b[t) in the following forms: 



1. a(t) = A [p 4- e at ) m with A> 0, a > 0, /3>0, m > 1 (emergent expansion) ]35l. 1431 ] 



2. a(t) = exp{A{hit) rn ) with Am > 0, m > 1 (logamediate expansion) 
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42]. 



3. b{i) = exp(Bt n ) with B > 0; < ra < 1 (intermediate expansion) [361. l42f| . 



4. b(t) = exp(B(lnt) n ) with Bn > 0, n > 1 (logamediate expansion) ]36l. 142 ] . 

5. b(t) oc t n with n > 1 (power law expansion) [2( . 

Some authors first choose the scale factor in power law, exponential or in other forms and 
then find out other variables with some conditions under these solutions. This 'reverse way of 
investigations had earlier been used extensively by [37J] who chose various forms of scale factor and 
then found out the other variables from the field equations. Subsequently, this approach has been 
adopted by [380, who clearly stated "This is not the ideal way to find out the dynamics of the 
universe, as here the dynamics is assumed and then the fields are found out without any reference 
to the origin of the field. But in the absence of more rigorous ways, this kind of investigations 
collectively might finally indicate towards the path where one really has to search". In another 
study, reference [3lJ assumed scale factor in the power law form to model the potential by an 
inverse square law in terms of the tachyon field. 
In the present work we have adopted this "reverse approach" to reconstruct the potential of the 



tachyonic field in the anisotropic universe. Subsequently, we have investigated how the statefinder 
parameters {r, s} behave in this case. We have chosen combinations of choices of scale factors 
and investigated different scenarios under the assumption that the universe is filled with tachyonic 
field. The basic equations are discussed in the following section. 



II. BASIC EQUATIONS 

In the present work we consider homogeneous and anisotropic iV-dimensional space-time model 
described by the line element [4 ll] 

d 
ds 2 = ds 2 FRW + Y J b2 ( t ) dx ^ ( 2 ) 

where d is the number of extra dimensions {d = N — 4) and ds 2 FRW represents the line element 
of the FRW metric in four dimensions is given by 



ds FRW = - dt2 + a2 W 



' /r> +r 2 (^ 2 +sin 2 ^ 2 ) 



(3) 



1 — kr 2 

where a(t) and b(t) are the functions of t alone represent the scale factors of 4-dimensional 
space time and extra <i-dimensions respectively. Here k (= 0, ± 1) is the curvature index of the 
corresponding 3-space, so that the above Universe is described as flat, closed and open respectively. 

The Einstein's field equations for the above non-vacuum higher dimensional space-time symme- 
try are 



,'d 2 + k\ D d 2 b 2 db 2 
3 (— )=D-Y¥ + 8¥ + P (4) 
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and 



i db Db h 2 p 

b + \b + Db-b 2= -2 (6) 
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Fig. 1 shows the variations of V against </>, for 

A = 1, B = 2, k = — 1, 0, 1, a = 2.5, /? = 3, to = 2, n = 3, d = 5 and Fig. 2 shows the variation of the 

statcfinder parameters r against s for A = 2, B = 3, a = 3, /3 = 2, m = 3, n = 2, d = 5 in the case of 

Tachyonic field for Emergent-Powerlaw Scenario. 



where p and p are energy density and isotropic pressure respectively. Here we choose here 
8irG = c = 1 and D 2 = b d (t), so we have % = f | and % = \ \ + ^f^g. 



Sahni et al [34|] proposed the trajectories in the {r,s} plane corresponding to different cosmo- 
logical models to depict qualitatively different behavior. Since this parameters are dimensionless 
so they allow us to characterize the properties of dark energy in a model independently. The 
statefinder diagnostic along with future SNAP observations may perhaps be used to discriminate 
between different dark energy models. The above statefinder parameters are given by [34[ | : 



H H , r-1 

1 + 3 7P + m and s = %3Ty 



(7) 



where q is the deceleration parameter defined by q = — 1 — jtj. In the anisotropic universe under 
consideration, the Hubble parameter H is given by H = ^3(3^ + d^) . In this paper we have 
reconstructed the statefinder parameters for different scenarios to be discussed in the subsequent 
sections. Finally we graphically analyzed {r, s} trajectory for various scenarios. 



Emergent-Power law Scenario: 




Fig.3 
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Fig. 3 shows the variations of V against <f>, for A = 2, B = 1, k = — 1, 0, 1, a = 3,y3 = 2, m = 3, n = .8, d = 3 

and Fig. 4 shows the variation of the statefinder parameters r against s for 
A = 2, B = 3, a = 2, (3 — 7,m — 2,n — 0.5, d = 5 in the case of Tachyonic field for Emergent-Intermediate 

Scenario. 



In this section we consider a combination of scale factors a and b in the emergent and power laws 
form respectively. As we are considering anisotropic universe, this combination appears feasible. 
We name this scenario as "emergent-power law scenario". In this scenario we consider a(t) and 
b(t) as follows: 



a(t) = A((3 + e at ) m and b(t) = Bt n 
Using equations (4)-(6), we can find the expressions for V((p) and (j) as 

V(t) 
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Fig.6 



Fig. 5 shows the variations of V against </>, for 

A = 2, B = 1, k = —1,0, l,a = 0.1, /3 = l,m = 0.3, n = 5,d = 3 and Fig. 6 shows the variation of the 

statefinder parameters r against s for A = 2, B = 3, a = 3, /3 = 2, m = 3.5, n = 2, <i = 5 in the case of 

Tachyonic field for Emergent-Logamediate Scenario. 
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Using equation (7), we can find the expressions for r and s as 



3(d + 3)(-3e at mt 2 a 2 /3 + dn(e a< + (3) 2 ) (3 + d) 2 (3e Q4 mi 3 a 3 /3(-e Qi + /3) + 2dn(e at + /3) 3 ) 
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(11) 
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Fig. 7 shows the variations of V against (/>, for A = 2, B = 1, k = — 1, 0, 1, m = 10, n = 0.9, d = 3 and Fig. 8 

shows the variation of the statefinder parameters r against s for A = 2, B = 3, m = 8, n = 0.8, d = 5 in the 

case of Tachyonic field for Logamediate-Intermediate Scenario. 



The above forms of statefinder parameters pertain to the emergent-power law scenario. The 
statefinder parameters derived above are plotted as r — s trajectory in the figure 2 and the 
potential V corresponding to this scenario is plotted against the scalar field eft in figure 1. 

• Emergent-Intermediate Scenario: 



In this section we consider a combination of scale factors a and b in the forms of emergent and 
intermediate respectively. We dub this scenario as "Emergent-Intermediate Scenario". Thus, to 
describe this scenario we consider a(t) and b(t) as follows: 



a(t) = A ((3 + e at ) m and b(t) = exp(Bt n ) 



(13) 



Using equations (4)-(6), we can find the expressions for V{(j>) and 4> as 



i / p -2A(\nt) m 

V{t) = -J ((24fct 2 - Bde^^^n^i-i + n(4 + B(d + l)t n )))(lnt) 2 + 24A 2 e 2A ^ t ) m m 2 (lnt) 2m ) 

O y v I 111 L j 



x \/(-8e- 2A ( lni ) m A;t 2 + B 2 (-l + d)dn 2 t 2n + 4Am(lnt) m - 2 (4 - 4m + (4 + Bdnt n ){lnt) - 6Am(lnt) m )) 

(14) 



e 2A(int)-n n t) 2 (-8e- 2A ( ln *) m M 2 + S 2 (-l + d)dn 2 i 2n + 4Am(lnt) m - 2 
0(t) = / W1 + 



(24/bt 2 - Bde 2A ( hlt ) m nt n (-4 + n(4 + B(d + l)t n )))(ln t) 2 + 24A 2 e 2A ( lni ) m m 2 (ln£) 
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x (4 - 4m + (4 + Bdnt n ) (In t) - 6,4m(lnt) m )) 

G>6 (10) 

Using equation (7), we can find the expressions for r and s as 

3(d + 3)(Bd(-l + n)nt"(lnt) 2 + 3Am(-l + m - ln£)(lnfp) 



r = 1 + 



+ 



(Bdnt n In i + 3Am(ln t) m ) 2 

(3 + d) 2 (Bdn(2 + (-3 + n)n)t n (ln t) 3 + 3Am(ln t) m (2 + (-3 + m)m + lni(3 - 3m + 2 In*))) 

(Bdnt n In t + 3Am(ln t) m ) 3 

(16) 

•((3 + d)(3(Bdnt n hit + 3Am(lnt) m )(5d(-l + n)nt n (lnt) 2 + 3Am(-l + m - lnt)(lnt) m ) 

QfRWn^W I qX^WV^wi , (d+3)(gd(-l+n)nf"(lnt) 2 +3Am(-l+m-lnt)(lnf) m ) x 
6{±fant mt-t6Am{mt) ) (^ H (Bdnt™ In t+3Am(lnt)™)2 J 

+(3 + d)(Bdn(2 + (-3 + n)n)t n (lnt) 3 + 3,4m(lnt) m (2 + (-3 + m)m + lnt(3 - 3m + 2 hit)))))" 

(17) 

The forms of potential and r — s trajectory derived in this scenario are presented in figures 3 and 
4 respectively. 

• Emergent-Logamediate Scenario: 

In this section we consider the scale factor a in emergent form and b in the logamediate form. 
We dub this scenario as "Emergent-Logamediate Scenario". Thus, in this scenario we consider ait) 
and b(t) as follows: 

a(t) = A (/? + e at ) m and b(t) = exp(B(lnt) n ) (18) 
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Using equations (4)- (6), we can find the expressions for V((p) and 4> as 
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Using equation (7), we can find the expressions for r and s as 
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In figures 5 and 6 we have plotted the tachyonic field potential against V against scalar field <f> for 
this scenario and r — s trajectory respectively. 

• Logamediate-Intermediate Scenario: 
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In this section we consider a in the logamediate form and b in the intermediate form. We name 
this scenario as "logamediate-intermediate Scenario". Hence, in this scenario we consider a(t) and 
b(t) as follows: 

a(t) = exp(A(}nt) m ) and b(t) = exp(Bt n ) (23) 

Using equations (4)- (6), we can find the expressions for V((p) and as 



v(t) = I 



— 2A(\n t) m 

x i I ~ 6 An xo ((24fct 2 - Bde 2A ( lnt ) m nt n (-A + n(4 + B(l + d)i n )))(lni) 2 + 24A 2 e 2A ( ln ') m m 2 (lni) 2 "^ 
t 4 (mt) z 

^J {-Se~ 2A ^ nt ) m kt 2 + B 2 (-l + d)dn 2 t 2n + 4Am(lni) m - 2 (4 - 4m + (4 + 5dnt n ) hit - 6,4m(lnt) m )) 

(24) 



, / e 2A(lnf)-( lnt )2(_ 8e -2A(lnt)"' A . t 2 + j g2(_ 1 + d ) dn 2 t 2n + ^m^t)^" 2 

° (/) ' A ' + (24/fct 2 - Bde 2 ^( lni ) m n^(-4 + n(4 + 5(1 + d)t n ))){\nt) 2 + 24A 2 e 2A ( lni ) mm2 ( ln *) 2m 



x (4 - Am + (4 + Bdnt n ) In i - 6Am(ln t) m ) 
Using equation (7), we can find the expressions for r and s as 

3(d + 3)(Bd(n - l)nt n (lnt) 2 + 3Am{-lnt + m - l)(lni) m ) 



r = 1 + 



(Sdnt n In i + 3,4m(lnt) m ) 2 

((d + 3) 2 (5dn(2 + (n - 3)n)t"(lnt) 3 + 3,4m(lnt) m (2 + (m - 3)m + Int(3 - 3m + 2 hit)))) 

(5dnt n In t + 3Am(ln t) m ) 3 

(26) 

_((d + 3)(3(Bdnt n Int + 3Am(lnt) m )(Bd(-l + n)nt n (lnt) 2 + 3Am(-lnt + m - l)(lnt) m ) + (d + 3) 
3(i?dntMnt + 3^m(lnt)-)3(| + 3 ( d+3 )^"- 1 ^ n ( lnf ) 2+3Am (- fat+m - 1 )( lnf ) m ^ 



-2 t (Bdnt n In t+3Am(ln t) m ) 

x(Bdn(2 + (-3 + n)n)t n (lnt) 3 + 3Am(lnt) m (2 + (-3 + m)m + lnt(3 - 3m + 2 lnt))))) 



(27) 



In figures 7 and 8 we have plotted the tachyonic field potential against V against scalar field 
for this "logamediate-intermediate scenario" and r — s trajectory respectively. 
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III. DISCUSSIONS 

In this work, we have considered N (= 4 + d)-dimensional Einstein's field equations in which 
4-dimensional space-time is described by a FRW metric and that of the extra <i-dimensions by 
an Euclidean metric. We have considered four scenarios, namely, Emergent-Powerlaw, Emergent- 
Intermediate, Emergent- Log amediate and Log amediate- Intermediate scenarios where the universe 
is filled with Tachyonic Field dark energy. The natures of the potentials as well as dynamics 
of scalar fields for the Tachyonic Field dark energy models have been analyzed. The statefinder 
parameters have been considered and their natures have been investigated for four reconstructed 
scenarios of the universe. 

In the case of Emergent-Power law scenario, we have reconstructed particular forms of scale 
factors a and b where a belongs to Emergent Scenario and b belongs to Powerlaw scenario in 
such a way that there is no singularity for evolution of the anisotropic Universe. We have found 
(ft and potential V in terms of cosmic time t for Tachyonic Field. Here we have shown that 
the Emergent-Powerlaw scenario is possible for open, closed or flat Universe if the Universe 
contains Tachyonic Field. From figure 1 it has been seen that the potential always increases with 
Tachyonic Field. The {r, s} diagram (fig. 2) shows that the evolution of the Emergent-Powerlaw 
Universe starts from asymptotic Einstein's static era (r — > oo, s — > — oo) and goes to ACDM 
model (r = 1, s = 0). It is also observed that r, s are independent of the dimension d. So, from 
statefinder parameters, the behavior of different stages of the evolution of the Emergent-Powerlaw 
Universe have been generated. 

In the case of Emergent- Intermediate scenario, we have reconstructed particular forms of scale 
factors a and b where a belongs to Emergent Scenario and b belongs to Intermediate scenario 
in such a way that there is no singularity for evolution of the anisotropic Universe. We have 
found (ft and potential V in terms of cosmic time t for Tachyonic Field. Here we have shown 
that the Emergent-Intermediate scenario is possible for open, closed or flat Universe if the 
Universe contains Tachyonic Field. From figure 3 it has been seen that the potential always 
increases with Tachyonic Field. The {r, s} diagram (fig. 4) shows that the evolution of the 
Emergent-Intermediate Universe starts from asymptotic Einstein's static era (r — > oo, s — > — oo) 
and goes to ACDM model (r = 1, s = 0). It is also observed that r,s are independent of the 
dimension d. So, from statefinder parameters, the behavior of different stages of the evolution of 
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the Emergent-Intermediate Universe have been generated. 

In the case of Emergent- Log amediate scenario, we have reconstructed particular forms of scale 
factors a and b where a belongs to Emergent Scenario and b belongs to Logamediate scenario 
in such a way that there is no singularity for evolution of the anisotropic Universe. We have 
found (ft and potential V in terms of cosmic time t for Tachyonic Field. Here we have shown 
that the Emergent-Intermediate scenario is possible for open, closed or fiat Universe if the 
Universe contains Tachyonic Field. From figure 5 it has been seen that the potential always 
increases with Tachyonic Field. The {r, s} diagram (fig. 6) shows that the evolution of the 
Emergent-Logamediate Universe starts from asymptotic Einstein's static era (r — > oo, s — > — oo) 
and goes to ACDM model (r = 1, s = 0). It is also observed that r,s are independent of the 
dimension d. So, from statefinder parameters, the behavior of different stages of the evolution of 
the Emergent-Logamediate Universe have been generated. 

In the case of Log amediate- Intermediate scenario, we have reconstructed a particular forms 
of scale factors a and b where a belongs to Logamediate Scenario and b belongs to Intermediate 
scenario in such a way that there is no singularity for evolution of the anisotropic Universe. 
We have found (ft and potential V in terms of cosmic time t for Tachyonic Field. Here we 
have shown that the Logamediate-Intermediate scenario is possible for open, closed or fiat 
Universe if the Universe contains Tachyonic Field. From figure 7 it has been seen that the 
potential always decreases with Tachyonic Field. The {r, s} diagram (fig. 8) shows that the 
evolution of the Logamediate-Intermediate Universe starts from asymptotic Einstein's static era 
(r —)■ oo, s —> — oo) and goes to ACDM model (r = 1, s = 0). It is also observed that r,s are 
independent of the dimension d. So, from statefinder parameters, the behavior of different stages 
of the evolution of the Logamediate-Intermediate Universe have been generated. 



IV. ACKNOWLEDGMENTS 

The authors are thankful to the reviewers for providing constructive and thoughtful suggestions 
to enhance the quality of the paper. The second author is thankful to the Inter-University Centre 
for Astronomy and Astrophysics (IUCAA), Pune for providing Visiting Associateship to him. 



14 



[1] Y-F. Cai et al, Physics Reports 493 (2010) 1. 

[2] E.J. Copeland, M. Sami, S. Tsujikawa, Internal. J. Modern Phys. D 15 (2006) 1753. arXiv:hep 
th/0603057. 

[3] T. Padmanabhan, Curr. Set. 88 (2005) 1057. 

[4] M. Li et al. larXiv:1103.R87(M [astro-ph.CO]. 

[5] B. Ratra and P. J. E. Peebles, Phys. Rev. D 37 (1988) 3406. 

[6] S. Nojiri, S. D. Odintsov, Gen. Relativ. Grav. 38(2006)1285. 

[7] V. Gorini, A. Kamenshchik, U. Moschella, Phys. Rev. D 67(2003)063509. 

[8] L. P. Chimento, Phys. Rev. D 69 (2004) 123517. 

[9] W. Zhao, Phys. Lett. B 655 (2007) 97. 

[10] N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. Lett. B 429 (1998) 263. 

[11] K. Poppenhaeger, S. Hossenfelder, S. Hofmann, M. Bleicher, Phys. Lett. B 582 (2004) 1. 

[12] T. Kaluza Preus. Acad. Wiss Fl 9669(1921). 

[13] O. Klein, A. Phys 37 895 (1926). 

[14] L. Randall and R. Sundrum, Phys. Rev. Lett 83 3370 (1999). 

[15] P. S. Wesson, World Scientific Space-Time-Matter, (1999). 

[16] M. R. Setare, Phys. Lett. B 653, (2007) 116. 

[17] M. R. Setare, J. Sadeghi, A. R. Amani, Phys. Lett. B 673 (2009) 241. 

[18] N. Arkani-Hamed, S. Dimopoulos, N. Kaloper and J.March- Russel, Nucl. Phys. B 567 189 (2000). 

[19] D. H. Lyth, Phys. Lett. B 448 191 (2002); J. M. Cline, Phys. Rev. D 61 023513 (2000). 

[20] A. Lukas, B. A. Ovrut and D. Waldram, Phys. Rev. D 61 023506 (2000). 

[21] H. A. Chamblin and H. S. Reall, Nucl. Phys. B 562 133 (1999). 

[22] D. Panigrahi and S. Chatterjee, Gen. Rel. Grav. 40 833 (2008). 

[23] J. Ibancz and E. Vardaguer, Astrophys. J. 306 401 (1986). 

[24] K. D. Krori, P. Borgohain and K. Das, Gen. Rev. Grav. 22 791 (1990). 

[25] R. J. Gleiser and M. C. Diaz, Phys. Rev. D 37 3761 (1988). 

[26] B. C. Paul, Phys. Rev. D 64 027302 (2001). 

[27] D. S. Gorbunov and S. M. Sibirvakov. rarXiv:0804.2248[ hep-th] . 

[28] G. Peng et al, Int. J. Theor. Phys. 41 1127 (2002). 

[29] D. Panigrahi and S. Chatterjee, Int. J. Mod. Phys. A 21 6491 (2006). 

[30] G. W. Gibbons, Phys. Lett. B 537(2002)1. 

[31] A. Feinstcin, Phys. Rev. D 66(2002)063511. 

[32] M. Sami, P. Chingangbam and T. Qurcshi, Phys. Rev. D 66 (2002) 043530; 



15 



M. Jamil, A. Sheykhi, Int.J. Theor. Phys 50 (2011) 625. 

[33] S. Chattopadhyay, U. Debnath, G. Chattopadhyay, Astrophys. Space Sci. 314 (2008) 41. 

[34] V. Sahni et al: JETP Lett. 77 (2003) 201. 

[35] S. Mukhcrjcc ct al. Class. Quantum Grav. 23 (2006) 6927. 

[36] J. D. Barrow, N. J. Nunes,Phys. Rev. D 76 (2007) 043501. 

[37] G. F. R. Ellis, M. Madscn, Class. Quantum Grav. 8 (1991) 667. 

[38] N. Banerjee, S. Das, Gen. Relativ. Grav. 37 (2005) 1695. 

[39] M. Sami, Mod. Phys. Lett. A 18(2003) 691. 

[40] Y. S. Piao, R. G. Cai, X. M. Zhang and Y. Z. Zhang, Phys. Rev. D 66(2002)121301. 

[41] S. Chakraborty and U. Debnath Int J Theor Phys 49(2010) 1693-1698. 

[42] S. Chakraborty and U. Debnath Int J Modern Phys A Vol. 25 24(2010) 4691-4701. 

[43] S. Chakraborty and U. Debnath Int J Theor Phys 50(2011) 80-87. 



